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Abstract 

Assuming so called global duality we argue that it is very likely that local 
duality needed to obtain results for the hadronic width of heavy meson de- 
cays within the 1/tmq expansion holds. Hence, if the discrepancy between 
experiment and the theory concerning charm counting, the semileptonic 
branching fraction and the lifetimes of b hadrons persist, it may be taken 
as a hint at some qualitatively new effect in (nonperturbative) QCD or 
even as a new physics. 



1 Introduction 



The 1 / rriQ expansion allows us to perform calculations of inclusive heavy hadron 
decay rates in a model independent and QCD based framework [|TJ. Unlike for 
exclusive decays we are even able to deal with purely hadronic processes using 
this method based on operator product expansion and Heavy Quark Effective 
Theory. This in turn allows us to calculate lifetimes and branching fraction for b 
hardrons within this QCD based approach. 

The leading order of the 1/w,q expansion is generically simply the parton 
model expression. The corresponding decay rates in the parton model are by now 
known for more than twenty years, and two problems with the rates calculated 
in this way have been noticed since then. Firstly, in the pure parton model, the 
semileptonic branching fraction comes out to be too large as compared to the 
measurements and secondly in the parton model all b hadron lifetimes are the 
same and equal the b quark lifetime. 

The hope was to fix these problems by including perturbative and nonper- 
turbative corrections originating from terms of order a s (mg) and l/m^l/mg. 
It turned out soon that the nonperturbative contributions are way too small to 
explain the semileptonic branching fraction || as well as the lifetime difference 
between the B mesons and the baryon |J. In fact the nonperturbative contri- 
butions are so small that they can safely be ignored given the present experimental 
situation. 

The purely perturbative corrections to the semileptonic decays have been 
calculated some time ago || and were also found to be too small to explain 
the data. The perturbative corrections to the hadronic decays have also been 
calculated in the late seventies for the case of vanishing masses of the final state 
quarks; again they are not sufficient to explain the semileptonic braching fraction 
(see e.g. refs. @, |2[ for the recent reviews). 

This was the motivation for Bagan and coworkers to perform the full 0(a s ) 
calculation including the mass effects of the final state c quarks in the channel 
b — > ccs. In fact, Bagan et al. find a substantial enhancement for this channel, 
which, together with the more modest enhancement in b — > cud channel due to 
the same mechanism almost explains the semileptonic branching fraction. This 
channel has also been discussed previously in this context in |7|], where it was 
argued that the operator product expansion could fail due to the small energy 
release in this channel, since the c quarks are heavy, or that parton hadron duality 
could fail due to resonances that are close to the point where the hadronic spectral 
function is replaced by the partonic one. 

However, one cannot blame the problem with the semileptonic branching frac- 
tion completely to the b — > ccs channel, since this would increase this channel by 
something like (30 - 40)%. This on the other hand contradicts charm counting in 
B decays, since this would give for the number n c of c and c quarks created per 
B decay n c ~ 1.3 which is about two standard deviations away from the experi- 
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mentally observed value of n c = 1.07 ± 0.07 ||. The QCD radiative corrections 
calculated in || also enhance this channel; the full 0(a s ) calculation also gives 
n c ~ 1.3, again in disagreement to charm counting. 

Another problem for the 1 / rriQ expansion is the lifetime ratio of the A;, baryon 
and the B meson. While the I/tuq expansion predicts a value in excess of 0.9 for 



the ratio r(A b )/r(-B + ), recent measurements yield r(A 6 )/r(5 + ) = 0.7 ±0.1 |L0 
which is a deviation at the level of two standard deviations. 

These two problems, namely the semileptonic branching fraction in combi- 
nation with charm counting and the lifetime of the A& baryon, have attracted 
considerable attention over the last time, although they are not yet of convincing 
significance. Neither the nonperturbative nor the perturbative corrections seem 
to be big enough to explain the problems 0. In particular, as it was mentioned 
above, although the recently completed calculation of the complete 0(a s ) correc- 
tion |5| is close to explain the discrepancy in the semileptonic branching fraction, 
it does it by the price of generating disagreement with charm counting. 

All the arguments used in the 1/rriQ expansion are based on parton-hadron 
duality and it has been suggested that a breakdown of duality is the explanation 
of these discrepancies. It has been pointed out that unlike in the calculation of 
the semileptonic branching fraction, where global duality is needed, in the calcu- 
lation of the hadronic width of the B meson local duality needs to be invoked: In 
the calculation of the semileptonic width only a weighted average of the sepectral 
function of hadronic operators is replaced by the corresponding average of the 
parton model expression, while for the hadronic width one replaces the hadronic 
spectral function at a fixed kinematic point (i.e. locally) by the partonic expres- 
sion. 

In this short note we shall assume that global duality holds and investigate to 
what extend the hadronic spectral function may deviate from the paron model 
expressions. Our arguments are necessarily qualitative, since it involves the non- 
perturbative sector of QCD. Intuitively one expects that the final state quarks 
my form resonance states which may lead to a hadronic spectral function differ- 
ent from the partonic one. We shall use simple resonance models to study the 
possible effects. 



2 Parton Hadron Duality 

We shall start from the effective Hamiltonian describing nonleptonic B decays, 
neglecting b — > u transitions and penguins 

H e ff = %V d> Y. [C 2 (M1-75)c)(C/7m(1-75)A/) (1) 

U=u,c 
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where D is the Cabbibo-rotated down type quark field 



V ud d + V* s s 



(2) 



d c = v; d d + v; sS 



and Ci and C 2 are the Wilson coefficients in the leading log approximation 



Ci 
C 2 



'a s {M ] 



w) 



7+ 



'a s {M ] 



w , 



' a s {M w ) V + ( a s (M w ) 
, a s (m 6 ) y 



(3) 



Here 
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7+ 



a s (m 6 ) 



'2 7 " 



33 - 2AT 7 

Using this Hamiltonian we may define the forward matrix element 
T( q \vq) = J d 4 xe-^(B(v)\H eff (x)H eff (0)\B(v)) 



(4) 



where v — Pb/ttib is the velocity of the B meson. 

We shall use the correlator T to study parton hadron duality. The analytical 
structure of T for complex vq and fixed q 2 is given by the intermediate states \X) 
that may be excited by i? e // 



T{q\vq) = J d^x {2Txf5\m B v~q-P x )\(B{v)\H eff \X)\ 2 . 



(5) 



The support of the function T may easily be calculated. For a fixed q 2 , T is 
nonvanishing in the interval 



vq < 



2m 



(q 2 + m 2 B — fi 2 ) and vq > — - — (q 2 + 2m| + /i 2 ) 



2m i 



(6) 



where \x is the mass of the lightest hadronic state which may appear in the sum. 
In our case this state is a Dn state where both the D and the ir are at rest, so 
fj, = m D + m w . 

In the parton model these intermediate states are three quark states at the 
tree level, and the function T may be calculated to be 



T(q 2 ,vq) 



'G 2 F mf 



[B(v)\bb\B(v))N c [C 2 + Cl + 



v 96tt 3 J 
V: d F 2 (m 2 JQ 2 ) + V; S F 3 (™ 2 JQ 2 ) 



2dC 2 

N r 



do 



(7) 
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where Q = rnsv — q. Furthermore, we have neglected CKM suppressed terms 
proportional to V us and V c a and 

F 2 (x) = (1 -x 2 )(l -8x + x 2 ) - 12x 2 ln(x)6(l - x) (8) 

F 3 (x) = v(l - Ux - 2x 2 - 12x 3 ) + 24x 2 (l - x 2 ) In —^9(1 - 4x), (9) 

1 — v 



v = y/l-Ax (10) 
Note that the functions F 2 and F 3 have support in the intervals 

vq < -^—(q 2 + m 2 B -m 2 c ). (11) 
Arris 



and 



vq< -^—(q 2 + m 2 B -4m 2 c ) (12) 
ZrriB 



respectively. 

Parton hadron duality means that the correlator T parton calculated in the 
parton model is "on the average" equal to the true hadronic correlator T hadron 

d(vq)T parton (q 2 ,vq)w(vq) = J d(vq)T hadron (q 2 ,vq)w(vq), q 2 fixed (13) 

where w is some smooth weighting function having support over some averaging 
interval. Close to the regions where only the Dit or J/\& K states contribute the 
partonic and the hadronic results have in general a quite different shape, and the 
interval over which one needs to average in ( |T3"D is larger than if we consider the 
two functions at higher invariant masses. THis situation has been analyzed for 



the case of the total cross section e + e~ — > hadrons in ITT . 

For the calculation of the total hadronic width we have to consider the corre- 
lator at the point where the invariant mass of the final state hadrons is equal to 
the B meson mass 

T(B -> hadrons) oc T(0, 0), (14) 

and the common folklore is that at this high invariant mass one may savely shrink 
the interval of averaging in ( |13D to a point and simply equate the partonic with the 
hadronic correlator. In other words, we shall assume that "global" duality in the 
sense of a large averaging interval holds and discuss under which circumstances we 
may shrink the averaging interval to a single point, thereby going from "global" 
to "local" duality. 

The situation in the nonleptonic decays has to be compared with the one in 
semileptonic decays |12|, |13j, where the hadronic contribution is the product of 
the two hadronic currents. In this case q may be identified with the momentum 
transferred to the leptons. While the analytical structure is quite similar (the 
lowest state being a D meson at rest) the integration over the neutrino momen- 
tum introduces a smooth averaging such that parton hadron duality is guaranteed 
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in almost all phase space for the charged lepton, except for the endpoint, where 
the energy of the charged lepton becomes maximal. Only if we would consider 
doubly differential distributions such that the invariant mass of the final hadronic 
state would be fixed, a similar situation would occur as in the nonleptonic de- 
cays, namely that the interval over which we average shrinks to a point. In this 
sense one expects that the calculation of the lepton spectrum as well as of the 
total semileptonic rate is more reliable than the calculation of the total nonlep- 
tonic rate, since in the semileptonic case only global duality is needed, while the 
nonleptonic case rests on the assumption of global duality. 



3 Resonance Model 

Motivated by this reasoning we have considered the correlator for the nonlep- 
tonic decays in the context of very simple model. The idea is to identify possible 
contributions which could generate some structure of the hadronic result which 
could result in a difference between T parton (0,0) and T hadron (0, 0) hence a viola- 
tion of local duality, although global duality holds. As pointed out above, the 
partonic tree level result is a smooth function of Q 2 = m 2 + q 2 — 2mf } (vq), since 
it is obtained from a three particle phase space for the three final state quarks. 
One possibility one may consider is that two of these final state quarks form a 
resonance, corresponding e.g. to the semiinclusive final processes B — > Di s \-\- 
anything or B — > J/^/+ anything. In order to study the qualitative features 
of this idea, we consider as a simple model for such a situation the effective 
Hamiltonian 

H' eff = ^U(b lfl (l- l5 ) q )r (15) 

where M is the field of a vector like resonance of mass M made of a quark- 
antiquark pair, and is its coupling strength to the left handed color singlet 
b — > q current, including all the CKM matrix elements. A similar model has 



been introduced by Stech and Palmer ||I4| to analyze semiinclusive nonleptonic 
B decays. 

It is a simple exercise to calculate the corresponding correlator. The resonance 
contribution to T near threshold turns out to be 

T(q 2 ,vq) = ^-f 2 {B(v)\bb\B(v)) (16) 

^(Q2 _ (M _ m<3 )2)(Q2 _ (M + m9 )2) 0( g2 _ (M _ ^2) 

Q 2 -M 2 + m 2 g + J ^((Q 2 -m 2 q ) 2 -M*) 

Contributions of other resonances to the spectral function are smooth since these 
resonances lie away from threshold, and we do not write it explicitly. Generically 
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the phase space contributes the factor involving the square root, and in the limit 
in which M, Q 2 ^> m q the r.h.s. of (|l^) is a smooth function even at threshold 



T{q\vq) = ^fi(B(v)\bb\B(v)) (17) 

(Q 2 - M 2 ) 2 (Q 2 + M 2 )6(Q 2 - M 2 ) 

In our case m q is just the mass of the light strange quark. This means that 
a scenario of a resonance as modelled in ( |l5l) such that the thereshold is close to 
the mass of the B meson does not introduce any rapid variations in the hadronic 
spectral function. Furthermore, since we assume global duality in the sense of 
(P~5|) , one finds that the hadronic and partonic spectral functions cannot be dras- 



tically different, at least not in a scenario as moddelsed by the Hamiltonian ([151) . 
Hence we conclude that it is quite likely that the partonic correlator T is a good 
approximation to the hadronic, even for the case of a resonance as we introduced 
above fairly close to the mass of the B meson. 

Typically the correlators obtained from intermediate states with two or more 
partons or hadrons are smooth functions of the variable Q 2 = m 2 B + q 2 — 2rrib(vq) 
and any resonance formation involving at least a two particle intermediate state 
does not introduce rapid variations of the correlator. For all these states we 
may apply the same argument as above, leading to the statement that it is quite 
likely that local duality holds. In this framwork, the only way to introduce strong 
variations is if there is a relatively narrow resonance R c leading to a single particle 
intermediate state in the correlator T. 

In order not to get into conflict with charm counting, this resonance state 
has to appear in the channel b —>■ cud, thereby enhancing this channel by some 
amount, which then would explain the semileptonic branching fraction without 
enhancing the charm created in B decays. In this case, the correlator would take 
the form 

T(,>,) = 2Mr: (Q2 _ M2)2 + M2r2 

where M and V are the mass and the width of R c respectively. If this resonance 
would be sufficiently narrow and its mass sufficiently close to the mass of the B 
meson, parton hadron duality could still hold on the average, but the interval of 
averaging would simply be too large to justify the calculation of the total hadronic 
width using local duality. 

However, this possibility is quite unlikely. First of all, it would not help to ex- 
plain the lifetime problem with the A& baryon, since such a resonance state would 
increase the width of the B mesons due to an increase in the channel b — > cud, 
and hence would make their lifetime shorter compared to the parton model value. 
Assuming that there is not such a resonance in the channel involving the baryons, 
this would increase the ratio r(A;,)/r(i?), which will make the discrepancy even 
stronger. 



.2 ^ _ ^^(B(v)\H eff \R c )(R c \H eff \B(v)} 
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Furthermore, the effect of such a resonance would be visible as well in the 
semileptonic channel, if the resonance state would have isospin / = 1/2. However, 
one may avoid any influence on the semileptonic decays, if the isospin of R c is 
I = 3/2. This would mean that in a quark model picture this state is a four 
quark state consisting of a two quarks and two antiquarks. 

However, the b — ► cud piece of the effective Hamiltonian ([!]) is an isospin 
triplet with I 3 = — 1, and hence the hadronic decays of the charged B mesons 
will have only final states with isospin / = 3/2, while the final states of the 
nonleptonic decays of the neutral B meson will be a mixture of I = 3/2 and 
I = 1/2. Consequently an enhancement of the I = 3/2 contribution through 
a resonance would make the lifetime of the charged B meson shorter than the 
one of the neutral states, which is in conflict with the data on the lifetime ratio 
between charged and neutral B mesons, which indicate that the lifetimes are 
equal t(B + )/t(B°) = 1.009 ±0.069 Hence there is not much room for such a 
resonance state and even if it would exist, its contribution would be not sufficient 
to explain the semileptonic branching fraction problem. 

4 Conclusions 

Even after more detailed calculations of the QCD radiative corrections the situa- 
tion concerning charm counting and the semileptonic branching fraction remains 
inconclusive. A failure of local parton hadron duality as it is needed to obtain 
results within the 1/rriQ expansion has been considered as a possible solution 
of the problem, but we have argued based on models that local parton hadron 
duality likely to hold, if global duality is assumed. 

Global duality, on the other hand, is a concept without which we could not 
relate many theoretical results to data, which in turn are compatible with the 
assumption of global duality. Hence the semileptonic branching fraction problem 
and the charm counting remain open questions in view of the present data. Ei- 
ther the data on charm counting change, such that n c ~ 1.3. In this case the 
problem is solved, and such a result would also be in agreement with the recent 
calculation of the full 0(a s ) corrections |J. However, if data remain as they are 
and the errors become smaller, than there have to be either yet unknown effects 
in nonperturbative QCD, or the discrepancies find an explanation in terms of 
physics beyond the standard model [16|]. However, we find an explanation 
in terms of QCD effects more likely, since physics beyond the standard model 
involves scales much larger than the B meson mass and hence will affect the 
decays of the b quark rather than the ones of the hadrons, in other words, such 
an explanation cannot be the source of the Af, lifetime problem. 
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